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Abstract 

Particle creation in spacetimes with a warped extra dimension is studied. In particular, we 
investigate the dynamics of a conformally coupled, massless scalar field in a five dimensional warped 
geometry where the induced metric on the 3-branes is that of a spatially flat cosmological model. 
We look at situations where the scale of the extra dimension is assumed (i) to be time independent 
or (ii) to have specific functional forms for time dependence. The warp factor is chosen to be that 
of the Randall-Sundrum model. With particular choices for the functional form of the scale factor 
(and also the function characterising the time evolution of the extra dimension) we obtain the 
the particle number and energy densities after solving (wherever possible, analytically but, 
otherwise, numerically) the conformal scalar field equations. The behaviour of these quantities 
for the massless and massive Kaluza-Klein modes are examined. Our results show the effect of a 
warped extra dimension on particle creation and illustrate how the nature of particle production on 
the brane depends on the nature of warping, type of cosmological evolution as well as the temporal 
evolution of the extra dimension. 
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I. INTRODUCTION 

BeEfinninEf with Kaluza-Klein hi, a large variety of models with extra dimensions have 

[in 

been proposed over the years. The recent brane- world models |2|-|5[ where our world is 
viewed as a four dimensional sub-manifold (3-brane) embedded in five dimensions are ac- 
tively pursued today largely because of their potential in proposing achievable experimental 
signatures of extra dimensions. Among braneworld models, the warped type necessarily 
assumes a curved higher dimensional spacetime and the line element on the 3-brane is 
scaled by a warp factor, thereby rendering the higher dimensional metric non-factorisable. 
The brane-world models seem to provide a viable resolution of the long-standing hierarchy 
problem in high energy physics. The warped braneworlds also suggest a dynamic way of 
compactification by proposing the idea of localisation of fields on the brane jol (and references 
therein). 

A study of quantum field theory in the context of warped spacetimes where there is an 
extra dimension is therefore not inappropriate. However, not much has been done along 
this direction. Leaving aside warping, the analysis of quantum field theory in a higher di- 
mensional spacetime with a Kaluza-Klein-like extra dimension has been discussed by some- 
notable among them being the works reported in [3], [sl, j^, [l^. More recently, Saharian 
fll| (and references therein) has discussed quantum fields in such spacetimes quite exten- 
sively (though, somewhat on the formal side) in a series of papers. Gravitational particle 
production in braneworld cosmology and its implications have been studied in 12|. 

The questions we address in this article are the following. Particle creation in a cosmo- 
logical background spacetime using the formalism of quantum fields in curved spacetime [l^ 
is a well-studied subject. How do the presence of extra dimensions affect particle creation? 
Further, how does a warping of the higher dimensional spacetime create differences, if any. 
The simple answer to the first question is related to the fact that now we do not refer to 
\Pkf but we must consider \l3k,k5f, where is the momentum associated with the extra 
dimension. The eigenvalues would have to be obtained by solving the equation for the 
extra dimensional part of the field (assuming separability) and the corresponding equation 
(say, for a scalar or a vector or a fermion). In the case of a KK like extra dimension k^ = ^ 
where is related to the radius of the compact extra space (say a circle). Not so when 
we have a warped spacetime. Here, the equation for the extra dimensional piece of the 
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field would be different for different types of warping and hence, the resulting solutions and 
eigenvalues will, obviously differ. Further, in a two brane scenario, appropriate boundary 
conditions need to be imposed and thus would take on only those allowed values, such 
that the boundary conditions are obeyed. The question therefore comes up: how do the 
values for k^^ as well as the nature of evolution of the scale factor and the extra dimensions 
affect particle creation? We shall provide illustrations to this end in the rest of this article. 
In addition, we also consider the case of a time-dependent extra dimension. In this con- 
text, we investigate how the nature of time evolution of the extra dimension affects particle 
creation characteristics. 

The plan of the article is as follows. In Section II, we discuss the conformally coupled 
scalar field equation and also the numerical method of solving the equations. The analytic 
formalism (for four dimensional cosmological spacetimes) is well-known and given in 13 1. 
Section III contains the analysis of the extra dimensional part (assuming separability) of 
the scalar field equation. Then, in Section IV, we find out the the number and energy 
density of the created particles for the four dimensional scenario with a massive, confor- 
mally coupled scalar field. Sections V and VI deal with the similar analysis for a massless 
scalar field in the presence of a time-independent extra dimension and a time dependent 
extra dimension, respectively. Section VII discusses the thermal/non-thermal nature of the 
spectra. Section VIII analyses zero mode particle creation and, finally, in Section IX, we 
conclude with comments and suggestions on future work. 



II. QUANTUM FIELD COUPLED TO A SPACETIME WITH AN EXTRA DI- 
MENSION 

Let us consider the background line element (using conformal time) to be generically of 
the form: 

ds" = e^f^''\^{r])[-drf + dx^ + V + dz''] + (j)''{v)da\ (2.1) 

where e^^^"'^ is the warp factor, a(r]) and (f){r]) are the scale factors associated with the 
ordinary space (x) and the extra dimension (a) respectively, rj denotes conformal time. 
A scalar field tlj{ri,x,a) conformally coupled to the above metric satisfies the following 
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Klein-Gordon equation, 



g-2/(a) Q2^ g-2/{a) /Q2^ Q2^ Q2^\ ^ Q2^ 



Sr]^ \ 9?/^ 



+ 



+ 



\ a (p I drj ' (j)^ da 



- (m^ + ^i?)^ = 0, 



(2.2) 



where m is the mass of the scalar particle, ^ is the conformal coupling constant and R is 
the five-dimensional curvature scalar. A dot ( ) denotes difi^erentiation w.r.t r] and a prime 
(/) denotes differentiation w.r.t a. 

We now concentrate on conformally coupled (^ = ^ in five dimensions) massless particles 
(m = 0). The scale factor evolution, the time-dependent extra dimension (we also discuss the 
time-independent case later) and the warp factor lead to distinct characteristics of particle 
creation. 

We separate variables using the following ansatz for the scalar field: 

1 



V^(r/,x,a) 



-Xi{v)F{^)G{a). 



The normalization condition for ip gives the Wronskian relation. 



(2.3) 



Xi Xi - XiXi = ^■ 



(2.4) 



Let 



and, 



F(x) [ dx^ 



G{a) ' G{a 
The above two assumptions imply. 



xiiv) ^ 

One can write Eg. (12.71) . as 
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Xiiv) = 0. 



(2.5) 



(2.6) 



(2.7) 



Uv)+ ^i+Q xi{v) = o 



(2.8) 



where, 



(2.9) 
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and 
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This equation admits WKB solutions of the form, 



Qldr) 



with a further restriction, 



where «/ and Pi are Bogohubov coefficients. 

Putting Eqs.([2H]) and fl2:T2D in Eqs.([2l]) and the condition ([23D, we get. 



and. 



i««r-iAr = i- 



?7— >oo 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



With the initial conditions, 

ai{r]o) = 1 and pi{r]o) = 0, 
the number of particles created in mode / (/ signifies both k and ka) is given by, 

Ni = lim lAI'. (2.16) 



Following Zel'dovich and Starobinsky [l^, Eqs. (l2.13l) and (12.141) . with the initial condi- 
tions (12.151) . can be cast in the form. 



dsi 




drj 




dvi 




drj 




dri 




drj 





1 + 2si) 



Ml 



Q 



2^1 



(2.17) 



with initial conditions, 



siivo) = ri{r]o) = vi{r]o) = 0, 



(2.18) 
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where, 



si = lA 



vi = 2Re \aiPle-^'r^"iv 



and 



ri = 2ilm laiPle-'^'T^i'^v . 



To get the number of particles created in mode / one has to solve the first order differential 
system f l2.17p with initial conditions (12.181) and determine si when t] oo. These equations 
can be evolved numerically using standard, easily available codes, in cases where one is 
unable to find an analytic solution. 

The number of created particles per unit volume, in ordinary space, in mode kcr is given 
at late times by 



III. THE ALLOWED VALUES OF k„ 

We first note the fact that there are specific allowed values of /Cq- which depend on the 
nature of the warp factor and the boundary conditions. To obtain these allowed values of 
kcr, we solved the Eq. (l2.6p . which is an eigenvalue equation for G(o"), for a typical functional 
form of f{a) (the RS solution). Other choices of f{a) can also be studied in a similar way. 

In Eq. (12.61) let us take G{a) = Gi{a)G2{(y), the requirement of a vanishing coefficient in 
the term involving G"^(cr) leads to the choice. 




(2.19) 



The energy density is given by 




(2.20) 



G2{a) = const, e'^^"). 



(3.1) 



and, subsequently, we have 
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2 
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+ A;^e-2^M Gi(a)=0. 



(3.2) 
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Now we can treat Eg. (13.20 as an eigenvalue equation for Gi{a) and investigate its solutions 
and the allowed values of k„. 

Following RS, we choose, /(a) = — 6 \a\. The Eq.f l3.2p . for a two-brane model, now 



becomes 

G'lia) + b {-6{a) + 6{a - r.vr)) - j + e^'^^\ 
The complete solution of this differential equation is. 



Gi{(y) = 0. (3.3) 



1 2b e-^'I'^l cos 



Gi{a) = Cd — n ,1.x ' where C is an arbitrary constant. (3.4) 

V cosyka-/ b) 



To crosscheck this result we use the following transformations, z = sgn{a) and G{z 



Gi{a) e '2' (due to [^]) in Eq. (13.31) . which leads to, 

+ kl G{z) = 0, (3.5) 

and confirms the previous result. Now, for a two-brane model, the boundary conditions at 
a = and cr = r^Tr imply 

dG{z) 3 b 



dz 2 1 + 6|;z|^^^^'-=o and 



(3.6) 



Eq. (13.51) with the above boundary conditions completely determine the allowed values for 
kfj through the following transcendental equation 



e + tan-i f --^e-^'^" 



tan _ , 

3 

We can approximate the above equation as 



-\6 where k„ = d b e-^'^^ . (3.7) 



2 

iaji9 = --e (3.8) 

for moderate values of br^ (i.e. ignoring the second term in the square brackets in the original 
transcendental equation). 

Solving graphically (see Fig. 1) we get ^ = 0, ±2.17, ±5.01, ±8.04, ±11.13, ±14.24, ±17.36 
and so on. We consider k^j = 9 in TeV units (as is normally done in warped braneworld 
scenarios) by absorbing the additional factor in the definition of the energy unit. 

In the case of an infinite extra dimension there may or may not be any discrete values of 
the modes. Even, if there are some (bound states within the volcano-shaped potential), it is 
clear that for a decaying warp factor of the RS type we would get a finite number of discrete 
values, and, beyond these values, we will have a continuum. In our calculations henceforth, 
we shall consider only the first three discrete values of k^ given above, for the two-brane 
model. 
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FIG. 1: 

IV. PARTICLE CREATION IN A 4D UNIVERSE 

Let us first look at the known analysis of a conformally coupled, massive scalar field in a 
four dimensional universe. In order to arrive at concrete results on quantities characterising 
particle creation we have to make a choice for the scale factor a{r]). This is chosen to be: 

a\v) = bl + bW (4.1) 

Note that the above choice gives a non-singular line element which has a bounce at rj = 
(FigH]). The approach to the big-bang singularity can be modeled using the limit bi —>■ 0. 
62 is known as the slowness parameter. In fact, for 61 — or 77 — ±00, the scale factor 
approaches that of a radiative universe. 

It is easy to show that in a four dimensional universe, the temporal part of a massive, 
conformally scalar field in the background given by: 

ds'^ = a{rif{-dr]'^ + dx^ + dy^ + dz^), (4.2) 

satisfies the equation 

Uv)+^Kv)Xi{v) = 0, (4.3) 

where 

n^{7]) = k^ + {bl + blr]^)m\ (4.4) 

As shown in [l^, one can construct an exact solution of the above equation in terms of 
parabolic cylinder functions. Using these solutions, one can easily obtain the particle number 
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density, 



exp 



—TT 



mh-: 



+ 



(4.5) 
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FIG. 2: Plots of and k^\(3k\'^ vs. k, respectively, in first and second columns, for the analytic 

solution in a universe without any extra dimension, for 3 values of m such as 2.17(long dashes), 
5.01(small dashes), 8.04(dotted). Here 62 = !> and in the 1st row, bi = \/0.2; in the 2nd row. 



61 = VOM and in the 3rd row, h = VOM. 

The main features of the particle number density (fc^|/5fcp) and energy density (A;^|/5fcp), 
can be derived from Figj2l are the following: 

(a) the k value, where or peaks, is independent of the parameter bi (for same 
m) but increases with increasing m. 

(b) The heights of all the peaks increases with decreasing bi. For higher values of 61, the 
lower modes dominate, whereas, when we lower bi (ie. we get closer to a metric with a 
singularity at rj = 0) the higher modes become more and more dominant. Also, it is known 
that massless particles (m = modes) will not be created, as the background universe is 
conformally flat (check this using m = in the formula for |/9fc|^). 

(c) The nature of the spectrum of created particles can be identified, following with 
that of a non-relativistic thermal gas of particles with momentum ^ at a chemical potential 



"I"a^ and temperature 2i^3-kB Boltzmann constant). 



V. PARTICLE CREATION WITH A TIME-INDEPENDENT EXTRA DIMEN- 
SION 

Having looked at the four dimensional scenario and also introduced the scale factor which 
we shall be working with throughout, we now move on to the case of a time-independent 
extra dimension. Note that this is different from the usual Kaluza-Klein scenario because of 
warping and also because of the choice of the extra dimensional space which could be finite 
or infinite. A fair amount of work on such brane cosmological models (the 3-braneworlds are 
now the so-called FRW branes) has been carried out since the inception of the braneworld 
idea 



0. 



A. Approximate analytic solution 

The analysis of a massless scalar field conformally coupled with a 5d metric (12. ip is now 
carried out with the choice ^(r/) = h\ (a constant equal to the minimum value of the scale 
factor a{ri)). The Eg. (12.70 takes the form: 



^2 _|_ "l + '^2 7 y2 \ , "l"2 



= 0, (5.1) 



bl 7 8{bl + bW) 

We have not succeeded in finding an analytic solution of the above equation. But, Eq. (l5.ip 
and Eq.( ]4.3p . are similar (modulo some redefinitions), if we consider that the particle concept 
has a meaning only at rj ^ ±oo. In this case, the second term in the square brackets in Eq. 
(5.1) can be ignored and an approximate solution can thus be found which is the same as 
that obtained in the scenario of a four dimensional universe. This leads to 



|/3fc|^ = exp 



fk%i k^bi 
\ Kb2 b2 



(5.2) 



Hence, the particle number density becomes dependent on the size of the extra dimension. 
This approximation implies the equivalence of m and i.e. one may imagine a massive 
field in 4d as a projection of a massless scalar field residing in the 5d bulk. The variations 
of total particle number density and total energy density w.r.t k for different parameter 
dependences are plotted in the left box of FigEl The following features can be noted: 
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(a) The k value where /c^|/3fcp or k^lPkl"^ peaks increases with decreasing bi for same kcr and 
it also increases with increasing k^r. 

(b) The heights of all the peaks increases with decreasing b^. 

(c) The higher modes dominate for lower bi values. 

(d) There will be no k„ = modes created (unlike the 4D case). Recall that the conformal 
invariance is not broken in the field equation since we have taken the Q{ri) factor to be 
neghgible. 



B. Numerical solution 



The case discussed in the previous subsection is now re-analysed using the Zel'dovich- 
Starobinsky equations given in the previous section (taking into account the extra factor 
Q{f])) and the corresponding variations are plotted in the right box of Figj3]in order to check 
how good the approximate analytic solution is. 
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FIG. 3: Left box shows plots of the analytic solution and the right box depicts numerically evaluated 
and A;^|/3fcp vs. k in two columns in a universe with a static extra dimension (with 62 = 1); 
for 3 values of k^, such as 2.17(long dashes), 5.01(small dashes), 8.03(dotted), where in 1st row, 
61 = 0.2; in 2nd row, 61 = 0.1 and in 3rd row, 61 = 0.05. 



The differences that arise between the results in this and those quoted in the previous 
subsection are: 
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(a) For same fco- and bi, peaks occur at different k values. 

(b) The presence of a static extra dimension decreases the amphtudes of total particle 
number density and energy density (this is a signature of the Q{ri) term). Here, it is obvious 
that zero mode particles will be created, but this phenomenon is addressed later. 

(c) The spectrum is no longer that of a thermal gas of non-relativistic particles. This is 
evident from the figure above where we show the k'^\Pkf and k^\Pkf for the approximate 
and the numerical analysis carried out in this section. Although it might seem from the 
plots of the number and energy densities that the spectrum is thermal, we note that this is 
not the case. The obvious reason behind this is hidden in the Q{ri) factor which breaks the 
conformal symmetry and thereby leads to a non-thermal spectrum. One might take this as 
a signature of the presence of the extra dimension. Our results with a time-dependent extra 
dimension which we shall turn to now, also exhibit the same fact. 

VI. PARTICLE CREATION WITH A TIME DEPENDENT EXTRA DIMENSION 



a(r]) 

— — 0ir],b3<l) 

(t>{r], b3>l) 

-0.4 -0.2 0.2 0.4 




FIG. 4: The above plot shows the variation of a{rj) and (/>(r/) (for two different choices of 63; we 
have taken 61 = 0.2, 62 = 1, 63 = 1.2 and 0.8 ) 

What happens when ^(r/) is time-dependent? Such extra dimensions have not really 
been studied much in the context of braneworld models. The unwarped scenario with time 
dependent extra dimension has been analysed in [17] with the motivation of obtaining an 
accelerating scale factor. Exact solutions of the Einstein equations with a non-constant 
0(77) and with physically motivated matter sources are rare and difficult to find, though 



some examples and their cosmological implications have been discussed in 



181 ]. However, in 



the Kaluza-Klein context time dependent extra dimensions which decay in time along with 
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the expansion of the universe have been dealt with in detail (see papers on Kaluza-Klein 
cosmology in fiol). 

Here, we assume the same form (as in the previous sections) for the scale factor a{ri). 
The (j){ri) is chosen to be a function of rj which approaches the constant value bi at the 
asymptotic infinities (±00). The form we choose is: 

'^'(^) = f'A?ff,2 - (6-1) 
FigH] shows the variations of the scale factor for typical parameter values. (f){r]) in both cases 
(63 > 1 and 63 < 1) goes over to 61 as ^ 00 (which was the size of the extra dimension 
considered in the previous section). 
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FIG. 5: Plots of numerically evaluated k'^\f3k\'^ and k^\/3k\'^ vs. k in two columns, in a universe with 
a dynamic extra dimension (with 63 = 1.2 and 62 = 1) for 3 values of A;^ such as 2.17(long dashes), 
5.01(small dashes), 8.04(dotted), where in 1st row, bi = 0.2; in 2nd row, 61 = 0.1 and in 3rd row, 
hi = 0.05. 

The particle number density and energy density are evaluated numerically (as before) in 
the presence of the above form of a time-varying extra dimension, for a typical value of 63 
(> 1) and plotted in FigjSl In this case, the size of the time-varying extra dimension always 
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remains less than the minimum size of ordinary space. Apart from the features pointed out 
earher, we can see that for 63 > 1, the contribution of the higher modes diminishes very 
quickly as compared to the previous cases. The reverse is seen to happen for 63 < 1, i.e. the 
higher modes tend to become more and more dominant with decreasing values of 63. Both 
these results are however dependent on the choice of bi which is taken to be large enough 
(away from the approach to a singular a{ri) (61 — > 0)). But, for a fixed bi (say bi = 0.2) 
the marked differences in the number and energy densities for 63 > 1 and 63 < 1 are clearly 
visible from the next plots. 

A. Comparison 

In order to figure out the similarities and differences arising out of a time-dependent 
and independent extra dimension, we now make a comparison of the results of the previous 
subsections. The differences are mentioned pointwise below. 

In Figini as the caption suggests, we have placed the salient features of the integrands 
required to compute the total particle number density and the total energy density in the 
four different situations (universe without extra dimension, with static extra dimension, 
and with a dynamic extra dimension of two different kinds). We use the same value of 
the parameters bi (determining the minimum size of the ordinary scale factor) and 62 (the 
slowness parameter) and the same k^r values (which are the first three non-zero solutions 
found in Section III) in consecutive rows. 

(a) The presence of a static extra dimension decreases the amplitudes for each mode, but 
the higher modes become relatively more dominating (i.e. relative contribution of higher 
modes w.r.t lower modes increases) in comparison with the case of a 4d universe. 

(b) We have two different kinds of time- varying extra dimension (FigS]). If the size of the 
extra dimension always remains less than the value equal to the static case and asymptot- 
ically reaches the same value (see 0(?7, 63 > 1)) in FigJl]), contributions of higher modes are 
significantly diminished and heights of the peaks are also decreased compared to the other 
cases. However, for the extra dimension with 0(?7, 63 < 1), an exactly opposite effect is seen 
to happen. 

(c) In the last three rows, the A;- values where k'^\(3k\'^ or k^lPkl"^ peaks are the same and less 
than that in the 4d case. 
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FIG. 6: Comparison of plots of and vs. k in two columns with, h\ = 0.2 and 62 = 



in a universe - 1st row: without extra dimension and m 



-; 2nd row: with static extra dimension; 



3rd row: with dynamic extra dimension and 63 = 1.2; 4th row: with dynamic extra dimension and 
63 = 0.8; for 3 values of k„ such as 2.17(long dashes), 5.01 (small dashes), 8.04(dotted). 



(d) Previously, in models where the extra dimension was taken to be hidden and a decreas- 
ing function of time ?!, any nonzero mode, if excited, was found to soon dominate over 
the redshifted /c-modes because of ever-increasing blueshift (their physical frequency asymp- 
totically reaches a value equal to k„/bi, but in the braneworld scenario bi can be much 
larger compared to those Planck-sized extra dimensions) resulting in a breakdown of the 
background cosmological model due to the possibility of a large back reaction. This prob- 
lem may be resolved here for a dynamic extra dimension with 63 > 1, where the zero- mode 
is more dominant and higher mode contributions are quickly suppressed i.e production of 
[k, ka 7^ 0] particles can be controlled using suitable values of 63. 
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VII. NON THERMAL NATURE OF THE SPECTRUM 



In the cases of the massive, conformally coupled scalar field in four dimensions and 
the analytic treatment of the time independent extra dimension (assuming Q{i]) = 0) the 
spectrum of created particles can be identified with that of a non-relativistic thermal gas 
of particles. However, including the Q{ri) factor we find that the spectrum deviates from 
the thermal nature-we might call it nearly thermal though. The same features appear when 
we consider time-dependent extra dimensions. In this section, we outline these features by 
attempting to fit (by least-square fitting method) the numerical data using thermal as well 
as non-thermal profiles characterised by a set of parameters. The results are shown below. 
In the case of a static extra dimension we attempt to fit total particle number densities (top- 
left plot in right box of Fig|3]) with the following expressions for thermal and nonthermal 
profiles of the particle number densities. 



The first one is essentially the same thermal profile as Eg. (15.21) with two new (fitting) 
parameters [A and B) introduced. The second one is a nonthermal profile with another 
parameter (n) introduced to take care of the non-thermality. The best-fit plots are shown 
in FigIT] and FigJH 

(a) Deviations of the datapoints from thermal fits are very prominent for lower values of fco-, 
whereas for higher and higher values of k„ data points seem to converge on the thermal fits. 
These features are clearer when we use a non-thermal fit which seems to be a better fit in 
all the three cases, but the extent of non-thermality decreases with increasing (i.e., value 
of n converges towards n = 2), this is due to the fact that the factor Q/Q contributes lesser 
and lesser for larger and larger k^j. 

(b) The fall-off of for large k, are much faster than the thermal fits. Again this is 
confirmed through non-thermal fits where n is found to be greater than 2. 

(c) To figure out a single possible fitting formula for the non-thermal profile with a unique 
set of parameters one needs to do a rigorous statistical analysis. The immediately apparent 
features of the parameters in our fitting functions are the following. With increasing k„ we 
find: A increasing, B showing an oscillatory behaviour and n converging towards 2. These 




(7.1) 




(7.2) 
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patterns can be useful in determining the nature of non-thermality-an aspect which seems 
to be crucial in distinguishing between the presence and absence of extra dimensions. 
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FIG. 7: Best-fit plots using a thermal profile for datapoints of A;^|/3fcp vs. k (top-left plot in right 
box of Fig l3] where bi = 0.2,62 = 1, for 3 values of ko-, such as 2.17(upper curve), 5.01(middle 
curve) and 8.03(lower curve)), and the corresponding best-fit parameter values. 
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FIG. 8: Best-fit plots for a nonthermal profile for datapoints of k'^\Pk\'^ vs. k (top-left plot in right 
box of Fig l3] where 61 = 0.2,62 = for 3 values of k^, such as 2.17(upper curve), 5.01(middle 
curve) and 8.04(lower curve)), and the corresponding best-fitting parameter values. 
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VIII. THE ZERO MODE 

In it is pointed out that, in a radiative universe fco- = modes or the massless particles 
in the four dimensional picture will not be created if we have a static extra dimension because 
this situation is equivalent to the problem of massless particle creation in a conformally flat 
four dimensional cosmology. This accidental recovery of conformal invariance happens only 
because of the particular choice of a{ri) (~ rj, in radiative universe) with a static extra 
dimension. Even in the case of a time-varying extra dimension conformal invariance can be 
imposed by making Q{i]) = (the trivial case is, when a{ri) = (f){ri), which is not of any 
interest). Otherwise, the conformal invariance is always broken for a metric of type Eg. (12.1 p 
and the zero mode particles will indeed be produced with cosmological evolution. In FigJHl 
features of zero mode particle creation in presence of a static extra dimension (continuous 
line) is compared with the other two situations with different kinds of extra dimensions 
(63 > 1: smaller dashes and 63 < 1: longer dashes). 




2 4 6 8 10 12 14 ^ 2 4 6 8 ' 10 12 14 



k k 

FIG. 9: Features of zero mode particle creation in presence of - static extra dimension (continuous 
line), dynamic extra dimension with 63 = 1.2 (smaller dashes) and dynamic extra dimension with 
63 = 0.8 (longer dashes); when bi = 0.2 and 62 = 1- 

The contribution of zero mode particles, in the static case, lies between those of the two 
dynamic cases (this feature is similar for the other modes). For these particular choices of 
parameters, the zero modes are the most dominating modes. But the plots of the integrand 
in the total particle number density, asymptotically reaches a particular value as A; ^ 
(for the zero modes), which is not the case for k^r ^ modes. Moreover, nothing similar to 
these modes (as the m = modes) exist in a four dimensional universe (without any extra 
dimension), since such a geometry is conformally flat. 
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IX. CONCLUSIONS AND COMMENTS 



In this article, we have investigated the possible effects of the presence of a warped extra 
dimension (static or dynamic) on particle production in a 4d universe. Prom the numerical 
analysis represented via the figures, it is easy to point out the characteristic features. 

• The warp factor plays its part through the values of which are the allowed modes 
arising because of the warped extra dimension, ka- will have discrete values in the case of a 
two brane model, otherwise it will be a continuum. 

• The location and heights of the peaks in the plots for the particle number density and 
energy density (as a function of k) decreases in the presence of an extra dimension. This 
is a feature of the breaking of conformal invariance and is manifest through the presence of 
the Q{r]) factor. 

• For the scenario with a time-varying extra dimension we have looked at two types 
distinguished through the parameter 63 in the expression for (p{r]). We note an increase in 
the heights of the peaks for 0(7^,63 < 1) and a decrease for 0(?7, 63 > 1)). To understand 
further general features of the breaking of conformal invariance in the context of particle 
creation one needs to see what happens for different kinds of models (i.e choices of 0(77) and 
0(77) which yield different forms for Qijj)). 

• It may be mentioned that though these features look to be case specific and model 
dependent, it is worth noting that wc did not choose ^(r]) arbitrarily. Rather, as mentioned 
before, it is chosen to be equal to hi (avoiding the introduction of another new parameter) 
or asymptotically reaching the same value (in the dynamic cases). Moreover, this value is 
the minimum value of a{r]) during its evolution. In effect we are assuming some correlation 
between aijj) and ^(r/), i.e. there may be a single mechanism which drives both the scale 
factors and therefore their ought to be a connection between them. This aspect can become 
clearer if we try to find out such solutions of the Einstein equations with a driving source 
in the bulk. 

• We have also found a mechanism (for ^(77, 63 > 1)) to have control over contributions 
of higher non-zero k^ modes to the energy density of the universe. We can suppress the 
production of those higher modes at our will by choosing higher and higher value of 63. 
The back reaction thus remains negligible and the classical cosmological model which is our 
background for all these studies does not break down. 
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• To get an idea about the thermal/non-thermal nature of the spectrum of created 
particles we have tried to fit our numerical data with a thermal profile with a new set of 
parameters. Our analysis reveals that the spectrum could be assumed to be nearly thermal 
though marked deviations seem to appear for larger values of k. 

• As a next step, one should study particle creation in the context of more realistic 
cosmological models and also further investigate the possibility of creation of other types 
of particles (like fermions, or spin one particles etc.) on the brane. Also, studies on par- 
ticle production along the brane could provide insights into the recently proposed idea of 
braneworld isotropisation 20 1. 

We admit, in conclusion, that our results are based on a toy scenario and can only 
be thought of as a build-up towards the study of more relevant and realistic situations in 
future. However, we do believe that some of the features observed (eg. possible avoidance 
of significant back-reaction effects and non-thermal nature of the spectral profile) here will 
indeed be carried over in generically similar situations and can provide pointers towards a 
better understanding of quantum fields in the presence of a warped extra dimension. 
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